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ABSTRACT 

^ ■ We report the detection of CO ( J=3— >2) line emission from all three multiple 

images (A,B and C) of the intrinsically faint (~ 0.8 mJy) submillimeter-selected 
galaxy SMM J16359+6612. The brightest source of the submm continuum emis- 
sion (B) also corresponds to the brightest CO emission, which is centered at 
2=2.5168, consistent with the pre-existing redshift derived from Ha . The ob- 
q ■ served CO flux in the A, B and C images is 1.2, 3.5 and 1.6 Jy km s" 1 respectively, 

with a linewidth of 500 ± 100 km s~ x . After correcting for the lensing amplifica- 
tion, the CO flux corresponds to a molecular gas mass of ~ 2 x 10 10 M Q , while 
Q H - the extent of the CO emission indicates that the dynamical mass of the system 

O ' ~9xlO lo M . Two velocity components are seen in the CO spectra; these could 

^ ■ arise from either a rotating compact ring or disk of gas, or merging substructure. 

The star formation rate in this galaxy was previously derived to be ~ 100-500 M 
yr _1 . If all the CO emission arises from the inner few kpc of the galaxy and the 
^ ■ galactic CO-to-H 2 conversion factor holds, then the gas consumption timescale is 

a relatively short 40 Myr, and so the submm emission from SMM J16359+6612 
may be produced by a powerful, but short-lived circumnuclear starburst event in 
an otherwise normal and representative high-redshift galaxy. 
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1. Introduction 

When did the bulk of stars in galaxies form? A huge number of studies have sought 
to answer this question (e.g., Smail et al. 2002; Giavalisco et al. 2004; Steidel et al. 2004 
and references therein). Observational evidence indicates that the star formation activity 
and AGN activity peaked at z >1, very likely between 2 < z < 3. Sub- millimeter (submm) 
surveys have shown that a large fraction of the total energy output of galaxies in this critical 
epoch of galaxy formation comes from luminous, dusty galaxies (Smail, Ivison & Blain 
1997; Blain et al. 1999). The most luminous submm sources are widely believed to be the 
progenitors of today's massive ellipticals (Lilly et al. 1999; Smail et al. 2002). The slope 
of the submm source counts rises steeply, indicating strong evolution, and the background 
radiation intensity at submm wavelengths is believed to be dominated by sources with fluxes 
~1 mJy (Blain et al. 2002). To understand galaxy evolution, it is therefore imperative to 
study the typically faint but numerous sub-mJy submm galaxies (SMGs). Unfortunately, 
even with the best submm continuum data, these intrinsically faint sources are difficult, if 
not impossible to study at present. The deepest existing 850 ^m continuum surveys are all 
limited by confusion at 2 mJy (Blain et al. 2002). 

Over 400 SMGs have been identified so far, yet their distribution and evolution as a 
function of redshift has been hampered by difficulties in identifying optical counterparts and 
obtaining reliable redshifts until relatively recently (Chapman et al. 2003, 2004). At the 
faint end of the luminosity function (with intrinsic fluxes < 10 mJy), only five redshifts 
for submillimeter galaxies have been confirmed absolutely by the detection of CO spectral 
lines (Frayer et al. 1998, 1999; Neri et al. 2003). 1 This paper presents the detection of CO 
(J=3^2) emission using the Owens Valley Radio Observatory's Millimeter Array (OVRO 
MMA) from the faintest submm galaxy to date, the remarkable object SMM J16359+6612 
(Kneib et al. 2004) that lies in a rare high-magnification configuration behind the interven- 
ing massive cluster gravitational lens Abell 2218. Its intrinsic 850/im continuum flux after 
correcting for the well-constrained gravitational lensing effect of the foreground cluster is 
only ~0.8 mJy. 

2. OBSERVATIONS 

We used OVRO MMA to search for CO ( J=3^2) line emission from SMM J16359+6612 
at the redshift (z=2.5165) derived from Ha using Keck-NIRSPEC (Kneib et al. 2004) be- 



A more substantial sample of SMGs has recently been compiled by Greve et al. (2004) 
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tween October 2003 and April 2004 in three array configurations, L, E and H. The total 
integration time on source was 52.9 hours with a typical single-sideband temperature of 
250K. The digital spectrometer was configured to observe at 98.3297 GHz with a spectral 
bandwidth of 448 MHz and a spectral resolution of 4 MHz. The phase and pointing cen- 
ter was a(J2000) = 16:35:52.50 and 5(J2000) = +66:12:15.00. We used the nearby quasar 
J 1638+573 to correct for atmospheric phase and gain variations. The spectral bandpass was 
determined from observations of the quasars 3C84, 3C345 and/or 3c273. Uranus and/or 
3c273 were used for absolute flux calibration, leading to a 15% uncertainty in the overall 
flux scale. Data from each track was calibrated separately using the MMA software package 
(Scoville et al. 1993), and the «v-data were imaged using standard routines in the MIRIAD 
package (Sault, Teuben & Wright 1995). The resulting channel maps were binned to a spec- 
tral resolution of 16 MHz (48.8 km s _1 ), and have a synthesized beam of 4'.'53 x 4'/06 at a 
position angle of —59° with a noise level of 1.7 mJy beam -1 . The channel maps were then 
corrected for primary beam attenuation and combined into a velocity-integrated map (see 
Figure 1). In Figure 2 we show the CO spectra corresponding to each image, along with a 
typical representation of the noise. 



3. RESULTS 

The CO (J=3^2) emission is brightest from the central image. Co-adding the spec- 
tra from the peak emission in each image and fitting a single gaussian yields a redshift of 
z=2.5168±0.0003, consistent with the Ha redshift (z Ha =2.5165±0.0015; Kneib et al. 2004). 
The CO kinematics are intriguing. CO emission is present at > 2a over 400 km s _1 in two 
velocity components located approximately —160 and +120 km s _1 with velocity widths of 
AUfwhm ~ 200 km s" 1 each. The total CO linewidth AVf whm ~ 500 ± 100 km s -1 . A consis- 
tent profile can be seen in image C to the south. The northern image A shows only the red 
component of the CO emission but image A is at the edge of our primary beam and suffers 
from a lower signal to noise ratio. The CO flux, corrected for primary beam attenuation, 
in the three images is 1.2+0.1, 3.5+0.1 and 1.6+0.1 Jy km s -1 , respectively, for a total CO 
flux of 6.3 Jy±0.2 km s _1 . To measure the molecular gas mass in the system, we adopt 
a = 4 M (K km s _1 pc 2 ) -1 , similar to the Galactic value (see, e.g., Sanders, Scoville & 
Soifer 1991), WMAP cosmology (H =71 km s _1 Mpc -1 , fi M =0.27, fi vac =0.73), and lensing 
amplification numbers from Kneib et al. (2004). Our choice of a = 4 M (K km s _1 pc 2 ) -1 
is consistent with that estimated for Arp 220 (Scoville, Yun & Bryant 1997) after correct- 
ing for brightness ratio of T 6 [CO ( J=3^2)]/T 6 [CO (J=l->0)] ~ 0.6 typically observed in 
starbursts (Devereux et al. 1994). We follow the formulae in Solomon, Downes & Radford 
(1992) using a luminosity distance calculated using WMAP cosmology. We correct for the 



-4- 



effect of gravitational lensing by dividing the observed CO fluxes of each image by its lensing 
amplification, and then compare the results inferred from all three images. Since A is at the 
edge of the primary beam and suffers from lower signal to noise, we also give the masses 
derived for each image separately (Table 1). The total H 2 mass is ~ 2xl0 10 M . Note that 
a is uncertain and depends on metallicity and excitation; it may be as low as a = 1 M Q (K 
km s _1 pc 2 ) -1 (e.g., Solomon, Downes, Radford, & Barrett 1997), in which case the total H 2 
mass is ~ 5 x 1O 9 M , comparable to the molecular gas mass in the Milky Way. 

We estimate a dust mass of 2 x 10 8 M Q from the 850/zm flux density, assuming a dust 
temperature of Td = 45 K. The inferred gas-to-dust ratio of 100 is consistent with metal-rich 
galaxies in the local universe, suggesting that the central regions of SMM J16359+6612 have 
already been significantly enriched with metals. 

The CO emission from image B may be partially resolved in one dimension. Using the 
imfit task in MIRIAD we fit the velocity-integrated emission and measure a deconvolved size 
of 3'.'0 x 174 at a position angle of -68°. According to the Kneib et al. (2004) model, the 
linear stretch due to lensing amplification along this position angle is about a factor of 2. 
The major axis of the CO emission thus has an intrinsic extent of about 12kpc, while the 
emission in the direction of the minor axis is unresolved. Note that the size estimates are 
upper limits given the low signal to noise of the map; it is possible that the CO emission is 
originating from a significantly smaller area. 

4. A Circumnuclear Starburst In Progress? 

SMM J 16359+66 12 has several remarkable properties. The galaxy morphology indicates 
an extremely red central core (7 in Kneib et al. 2004) bracketed by blue features (a,/3). 
Unlike many submm galaxies, which show Lya is present in emission (e.g., Chapman et 
al. 2003), SMM J16359+6612 has Lya in absorption, as in most Lyman-break galaxies 
(LBGs) (Shapley et al. 2003). Given the complex morphology, the Lya absorption, the 
bright submm and infrared fluxes and compact Ha emission, the star-formation activity in 
SMM J16359+6612 is likely to be highly dust enshrouded. The size of the red core in HST 
images suggests that the star formation has an extent of only a few kpc (Kneib et al. 2004). 
The Ha linewidth in this galaxy is only ~280 km s _1 , which is greater than for most LBGs 
(e.g., Pettini et al. 2001; Erb et al. 2003), but less than typically found for high-redshift 
QSOs (e.g., Ivison et al. 1998). Nevertheless, the star formation rate estimated from the far- 
infrared luminosity may be as high as 500 M yr _1 (Kneib et al. 2004), an order of magnitude 
larger than in most LBGs (e.g., Pettini et al. 2001; Erb et al. 2003). Is SMM J16359+6612 
at a later evolutionary stage than submm sources previously studied in CO? Is it massive? 
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What will be the fate of this galaxy and how long can it maintain its starburst /possible 
AGN activity? 

The molecular gas observations offer answers to some of these questions. The molecular 
mass inferred from the CO data provides a modest reservoir of ~2xl0 10 M for this rate 
of star formation activity. The extent of the CO emission indicates that the molecular gas 
is likely located in the central few kpc, similar to the size of the red core in optical images 
(Kneib et al. 2004). If all the gas is consumed at the estimated current star formation rate 
of 100-500 M Q yr _1 , then the activity should cease in a mere 40-200 Myr. The nature of 
the star formation is reminiscent of circumnuclear starbursts in local infrared ultraluminous 
galaxies where the timescale for bursts is typically on the order of ~100 Myr, although the 
star formation rates and central molecular concentration in the lower redshift galaxies are 
both less by about a factor of 10. The gas and star formation activity in local starbursts are 
often distributed in a ring-like or disk-like morphology (e.g., Downes & Solomon 1998). In 
this respect too SMM J16359+6612 could resemble a circumnuclear starburst: the two CO 
velocity components may indicate a gas ring or disk morphology. If the CO were centrally 
concentrated (and in the absence of merging) we would expect a single broad velocity peak 
in the CO emission. 

There is one troubling discrepancy between the CO and Ha observations. The CO 
linewidth is about 500 km s -1 whereas the Ha linewidth is 280 km s _1 . One explanation 
for this difference is that the two lines trace different parts of the galactic potential: the Ha 
emission is likely confined to only the central 1-2 kpc in the source plane (Kneib et al. 2004) 
whereas the CO emission may come from a more extended region. This distribution could 
be consistent with a ring morphology for the CO distribution if the central molecular gas has 
been consumed by the starburst activity. Another possibility is that we are only observing 
the Ha emission from the near side of the galaxy or through windows where there is less 
extinction from dust. As the CO emission is only barely resolved, higher spatial resolution 
observations (e.g., with IRAM, CARMA 2 and eventually ALMA 3 ) are needed to understand 
the relative distribution of the molecular gas and star formation activity. 



2 The Combined Array for Research in Millimeter- wave Astronomy, http:/ /www. mmarray.org 
3 Atacama Large Millimeter Array, http://alma.nrao.edu 
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5. Is This Submm Source Fated to Become a Massive Early Type Galaxy? 

Do all galaxies pass through a bright-submm phase? This is an outstanding question 
in the literature (e.g., Lilly et al. 1999; Smail et al. 2002). Certainly the brightest submm 
sources at near-IR wavelengths are believed to be analogs of today's local, massive, early type 
galaxies (Smail et al. 2002, 2004). But what about the fainter submm sources that contribute 
the majority of the submm background? Do later Hubble types and less massive galaxies 
also pass through a bright submm phase when building their bulges and perhaps when their 
stellar disks begin to form? If so, then the corresponding fainter submm galaxies should not 
be very massive (< 10 12 M ). With the CO data for SMM J16359+6612, comparable to that 
available for cB58 (Baker et al. 2004a), we can make an attempt to address this fundamental 
question. 

Assuming that the molecular gas is tracing the galaxy mass at a radius of 6 kpc, and 
using the CO linewidth of 500 km s" 1 , we estimate a dynamical mass of 9xl0 10 M Q . Note, 
however, that the mass is uncertain by at least a factor of two and depends on several 
important assumptions. The mass may be smaller depending on the true size of the CO 
emitting region which is, at best, only barely resolved in these data. The mass could also be 
higher depending on the inclination of any disk present. 

Not surprisingly, the dynamical mass in SMM J16359+6612 is at the low end of SMGs 
(see review by Genzel et al. 2004 and references therein), and is higher than the typical 
dynamical masses (~1O 1O M ) inferred for LBGs (Pettini et al. 2001; Erb et al. 2003; Baker 
et al. 2004a,b). But unlike previous SMGs, and especially the strongly lensed LBG cB58 for 
which CO spectroscopy is available (Baker et al. 2004a), SMM J16359+6612 lacks a large 
molecular gas reservoir by fraction. The relatively low gas fraction suggests that the bulk of 
the stellar mass in this system is already in place and so the stellar mass is unlikely to increase 
significantly without subsequent merging activity. So does this imply that a submm-bright 
phase is typical for all galaxies including the less massive galaxies in the field today? 

SMM J16359+6612 is part of a close galaxy group: two other galaxies (#384/#478 and 
#273), neither of which are detected in the submm, are within 130 kpc, and all three are 
separated by less than 100 km s _1 from each other in redshift. The submm luminosity of 
SMM J16359+6612 may be boosted by star-formation triggered by the interaction with these 
neighbors (Kneib et al. 2004). This situation is analogous to local galaxies where otherwise 
normal spirals undergo massive bursts of star formation and molecular gas inflow upon an 
interaction or merger (e.g., Scoville, Yun & Bryant 1997 and references therein). We cannot 
conclusively state that the faint end of the submm population is composed exclusively of low 
mass systems. However, as SMM J16359+6612 appears to be less massive than the brightest 
submm sources, it is clear that a fraction of faint submm galaxies could be systems with low 
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gas fractions and more modest total masses that could have elevated submm fluxes due to 
intense but short-lived starburst activity. 

6. Conclusions 

Our detection of molecular gas from SMM J16359+6612, the faintest submm-selected 
galaxy observed in CO emission to date, reveals some remarkable properties. SMM J16359+6612 
is significantly smaller in extent, and has a lower dynamical mass and a lower gas mass frac- 
tion than previous SMGs. These factors all suggest that SMM J16359+6612 is unlikely to 
develop into a massive elliptical galaxy without extensive subsequent merging. The dis- 
tribution and kinematics of the molecular gas and star formation activity indicate that a 
circumnuclear starburst is likely to be responsible for the submm emission, in analogy with 
nearby ultraluminous infrared galaxies. The luminous activity is likely induced by interac- 
tions with two neighboring galaxies. SMM J16359+6612 indicates that at least a fraction 
of the sub-mJy submm sources are likely going through a bright submm phase due to an 
episode of intense star formation activity that is found in a galaxy of relatively modest mass, 
and triggered by environmental factors. This type of interaction-induced activity in a mod- 
estly massive system could help to explain the very strong clustering signal seen in submm 
galaxies Blain et al. (2004). 
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CO {3-2) del B50jjm eminum CO (3-2) on optical 




Fig. 1. — Contours show OVRO observations of the CO (J=3—>2) velocity-integrated, pri- 
mary beam corrected, emission overlaid on both 850yum emission (left panel) and an optical 
HST image (right panel). The three images are labeled A,B, and C following the convention 
of Kneib et al. (2004). The CO contour levels are -2,2,3,4,5... x 0.32 Jy km s™ 1 . The peak 
of the 850yum submm continuum emission coincides with B at 17 mJy 
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Velocity Integrated C0(3 — 2) emission 
from J16359+6612 behind A2218 
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Fig. 2. — A velocity-integrated, primary-beam corrected, CO (J=3^2) emission map of 
SMM J16359+6612. The small plus sign between B and C indicates the OVRO phase and 
pointing center. Hanning-smoothed spectra for the three lensed images are shown on the 
right and a noise spectrum is shown in a panel on the bottom left. The dashed line shows the 
2a hanning smoothed noise level. The CO and optically determined redshifts are indicated 
with arrows. 



Table 1. OVRO CO(3 -2) Data for the three images of SMM J16359+6612 



ID 


a(J2000) 


5(J2000) 


Sco 


Amplification 


Intrinsic Sco 


M(H 2 ) 








(Jy km s _1 ) 




(Jy km s _1 ) 


(10 10 M ) 


A 


16 A 35 m 54.814 s 


66°12'37" 


1.2±0.14 a 


14 


0.09" 


1.3 


B 


16 A 35 m 54.152 s 


66°12'24" 


3.5±0.12 


22 


0.16 


2.3 


C 


16 A 35 m 50.848 s 


66°12'06" 


1.6±0.13 


9 


0.18 


2.6 



a Image A is at the edge of the primary beam and suffers from low signal-to-noise ratio, as 
noted in the text. 



